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Abstract 

In this paper we overview some results on the hydrogen atom 
in external static uniform magnetic fields. We focus on the case of 
very strong magnetic field, B ^> Bq = 2.3 • 10 9 Gauss, use various 
approximate models and, particularly, in the adiabatic approxima- 
tion have calculated exactly the integral defining the effective poten- 
tial. This potential appears to be finite at z = 0. Our consider- 
ation of the problem of highly magnetized atoms and molecules is 
motivated by the recently developed MagneGas technology by San- 
tilli flhttp:/ /www.magncgas.com ). The ground state electron charge 



distribution of the hydrogen atom in an intense magnetic field is of 
a toroidal form, in agreement with that studied by Santilli. This 
physical picture is at the foundation of the new chemical species of 
magnecules proposed by Santilli. 







1 Introduction 



Weak external static uniform magnetic field B causes anomalous Zeeman 
splitting of the energy levels of the hydrogen atom, with ignorably small 
effect on the charge distribution of the electron. In the case of a more intense 
magnetic field which is strong enough to cause decoupling of a spin-orbital 
interaction (in atoms), ehB/2mc > AEjji — 10~ 3 eV, i.e. B ~ 10 5 Gauss, a 
normal Zeeman effect is observed, again with ignorably small deformation of 
the electron orbits. 

In the case of weak external magnetic field B, one can ignore the quadra- 
tic term in the field B because its contribution is small in comparison with 
that of the other terms in Schrodinger equation so that the linear approxi- 
mation in the field B can be used. In such a linear approximation, the wave 
function of electron remains unperturbed, with the only effect being the well 
known Zeeman splitting of the energy levels of the hydrogen atom. In both 
Zeeman effects, the energy of interaction of electron with the magnetic field is 
assumed to be much smaller than the binding energy of the hydrogen atom, 
ehB/2mc <C me 4: /2h 2 = 13.6 eV, i.e. the intensity of the magnetic field is 
much smaller than some characteristic value, B ^ B = 2.4 • 10 9 Gauss = 
240000 Tesla (1 Tesla = 10 4 Gauss). Thus, the action of a weak magnetic 
field can be treated as a small perturbation of the hydrogen atom. 

In the case of very strong magnetic field, B 3> Bq, the quadratic term in 
the field B makes a great contribution and can not be ignored. Calculations 
show that a considerable deformation of the electron charge distribution in 
the hydrogen atom occurs. Namely, under the influence of a very strong 
external magnetic field a magnetic confinement takes place, i.e. in the plane 
perpendicular to the direction of magnetic field the electron dynamics is 
determined mainly by the action of the magnetic field, while the Coulomb 
interaction of the electron with the nucleus can be viewed as a small per- 
turbation. This adiabatic approximation allows one to separate variables in 
Schrodinger equation |TJ. At the same time, in the direction along the di- 
rection of the magnetic field the motion of electron is governed both by the 
magnetic field effects and the Coulomb interaction of the electron with the 
nucleus. 

In this paper we briefly review some results on the hydrogen atom in very 
strong external static uniform magnetic fields, focusing on the basic physical 
picture derived from the Schrodinger equation. Our consideration of the 
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problem of atoms and molecules exposed to strong external magnetic field 
is related to MagneGas technology and PlasmaArcFlow hadronic molecular 
reactors recently developed by Santilli 0. 

The highest intensities maintained macroscopically at large distances in 
modern magnet laboratories are of the order of 10 5 ...10 6 Gauss (~ 50 Tesla), 
i.e. much below -Bo = 2.4 -10 9 Gauss (~ 10 5 Tesla). An extremely intense ex- 
ternal magnetic field, B > B s = B /a 2 = 4.4 -10 13 Gauss, corresponds to the 
interaction energy of the order of mass of electron, mc 2 = 0.5 MeV; a = e 2 /he 
is the fine structure constant. In this case, despite the fact that the extremely 
strong magnetic field does not make vacuum unstable in respect to creation 
of electron-positron pairs, one should account for relativistic and quantum 
electrodynamics (QED) effects, and invoke Dirac or Bethe-Salpeter equation. 
Such intensities are not currently available in laboratories. However, these 
are of interest in astrophysics, for example in studying an atmosphere of 
neutron stars and white dwarfs which is characterized by B ~ 10 9 . . . 10 13 
Gauss. 

In this paper, we shall not consider intensities as high as Schwinger value 
Bs, and restrict our review by the intensities 2.4 ■ 10 10 < B < 2.4 • 10 11 
Gauss (10£>o • • • 100£>o), at which (nonrelativistic) Schrodinger equation can 
be used to a very good accuracy, and the adiabatic approximation can be 
made. 

Relativistic and QED effects (loop contributions), as well as the effects 
such as those related to the finite mass, size, and magnetic moment of the 
nucleus, and the finite electromagnetic radius of electron, reveal themselves 
even at low magnetic field intensities, and can be accounted for as very small 
perturbations. These are beyond the scope of the present paper, while being 
of much importance in the high precision studies, such as those on stringent 
tests of the Lamb shift. 

It should be noted that a locally high-intensity magnetic field may arise 
in the plasma as the result of nonlinear effects, which can lead to creation 
of stable self-confined structures having a nontrivial topology with knots [[| . 
Particularly, Faddeev and Niemi || recently argued that the static equilib- 
rium configurations within the plasma are topologically stable solitons, that 
describe knotted and linked fluxtubes of helical magnetic fields. In the region 
close to such fluxtubes, we suppose the magnetic field intensity may be as 
high as B Q . In view of this, study of the action of strong magnetic field and 
the fluxtubes of magnetic fields on atoms and molecules becomes of much 
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interest in theoretical and applicational plasmachemistry. 

Possible applications are in MagneGas technology, theoretical founda- 
tions and exciting implications of which are developed by Santilli 0. We 
refer the reader to Ref. for the recent studies, description, and applica- 
tions of MagneGas technology and PlasmaArcFlow reactors. Remarkable ex- 
perimental results, including Gas-Chromatographic Mass-Spectroscopic and 
InfraRed spectrum data of Santilli's magnegas, clearly indicates the presence 
of unconventional chemical species of magnecules, which have not been iden- 
tified as conventional molecules, and are supposed to be due to the specific 
bonding of a magnetic origin. The basic physical picture of the magnecules 
proposed by Santilli is related to the toroidal orbitals of the electrons of 
atoms exposed to a very strong magnetic field. The results of the present 
paper demonstrate this physical picture on the basis of the Schrodinger equa- 
tion. 

As the result of the action of very strong magnetic field, atoms attain 
great binding energy as compared to the case of zero magnetic field. Even 
at intermediate B ~ Bq, the binding energy of atoms greatly deviates from 
that of zero-field case, and even lower field intensities may essentially affect 
chemical properties of molecules of heavy atoms. This enables creation of 
various other bound states in molecules, clusters and bulk matter |Lj, ]2|, §j|. 

The paper by Lai [|J, who focused on very strong magnetic fields, B 3> Bo, 
motivated by the astrophysical applications, gives a good survey of the early 
and recent studies in this field, including those on the intermediate range, 
B ~ B , multi-electron atoms, and H 2 molecule. Much number of papers 
using variational/numerical and/or analytical approaches to the problem of 
light and heavy atoms, ions, and H2 molecule in strong magnetic field, have 
been published within the last six years (see, e.g., references in [f|]). However, 
highly magnetized molecules of heavy atoms have not been systematically 
investigated. One of the surprising implications is that for some diatomic 
molecules of heavy atoms, the molecular binding energy is predicted to be 
several times bigger than the ground state energy of individual atoms. || 

2 Landau levels of a single electron 

To estimate intensity of the magnetic field which causes a considerable de- 
formation of the ground state electron orbit of the hydrogen atom, one can 
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formally compare the Bohr radius of the hydrogen atom in the ground state, 
in zero external magnetic field, clq = h 2 /me 2 ~ 0.53 • cm = 1 a.u., with 
the radius of orbit of a single electron moving in the external static uniform 
magnetic field B. Below, we give a brief summary of the latter issue. 




Figure 1: Classical orbit of the electron in an external static uniform magnetic 
field B = (0, 0, B) pointed along the z axis. The electron experiences a 
circular motion in the x, y plane and a linear motion in the z direction (a 
helical curve). 

The mean radius of the orbital of a single electron moving in a static 
uniform magnetic field can be calculated exactly using Schrodinger equation, 
and is given by 



where we have denoted 




(2.1) 



(2.2) 



B is the intensity of the magnetic field pointed along the z axis, B = (0, 0, B), 
f = (r, if, z) in cylindrical coordinates, and n = 0,1,... is the principal 
quantum number. Thus, the radius of the orbit takes discrete set of values 



(2.1), and is referred to as Landau radius. This is in contrast to the well 



known classical motion of electron in the external magnetic field (a helical 
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curve is shown in Fig. [I]), with the radius of the orbit being of a continuous 
set of values. The classical approach evidently can not be applied to study 
dynamics of the electron at atomic distances. 

Corresponding energy levels E n of a single electron moving in the external 
magnetic field are referred to as Landau energy levels, 

En = Ei + El = hQ(n + \)+^, (2.3) 

where 

n = — (2.4) 

mc 

is so called cyclotron frequency, and hk z is a projection of the electron's 
momentum hk on the direction of the magnetic field, — oo < k z < oo, m is 
the mass of electron, and — e is the charge of electron. 

Landau energy levels E^ correspond to a discrete set of round orbits 
of electron which are projected to the transverse plane. The energy e\ 
corresponds to the free motion of electron in parallel to the magnetic field 
(continuous spectrum), with a conserved momentum hk z along the magnetic 
field. 

Regarding the above presented review of Landau's results, we remind 
that in the general case of uniform external magnetic field the coordinate 
and spin components of the total wave function of the electron can always 
be separated. 

The corresponding coordinate component of the total wave function of 
the electron, obtained as an exact solution of Schrodinger equation for a 
single electron moving in the external magnetic field with vector-potential 
chosen as A r = A z = 0, A v = rH/2), 

~ ^ + \d r + 1^ + dl - 7 V + 2ry^) ^ = (2.5) 
is of the following form |l]: 

e ilip e ik z z 

1pn )S ,kA r ^i Z ) = J^llns{.ir 2 ) (2.6) 

y Ztx v Li 

where I ns (p) is Laguerre function, 

UP) = ^=e-^y n - s)/2 Qr*(p); (2.7) 
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Q™~ s is Laguerre polynomial, L is normalization constant, / = 0, ±1, ±2, . . . 
is azimuthal quantum number, s = n — I is the radial quantum number, and 
we have denoted, for brevity, 

p = 7r 2 . (2.8) 

The wave function ( |2.6| ) depends on the angle (p as exp[ilip], and the energy 
( |2.3|) does not depend on the azimuthal quantum number I, so the system 
reveals a rotational symmetry in respect to the angle ip and all the Landau 
orbits are round. 

Also, the energy ( |2.3|) does not depend on the radial quantum number 
s. This means a degeneracy of the system in respect to the position of the 
center of the orbit, and corresponds to an "unfixed" position of the center. 
At n — s > 0, the origin of the coordinate system is inside the orbit while at 
n — s < it is outside the orbit. The mean distance d between the origin of 
coordinate system and the center of orbit takes a discrete set of values, and 
is related to the quantum number s as follows: 



Clearly, the energy ( |2.3Q does not depend on this distance because of the 
homogeneity of the magnetic field, so it does not depend on s. However, 
it should be noted that the whole Landau orbit can not move on the (r, (p) 
plane continuously along the variation of the coordinate r. Instead, it can 
only "jump" to a certain distance, which depends on the intensity of the 
magnetic field, in accord to Eq. ( |2.9|) . 

The spin components of the total wave function are trivially 

V<-i) = (?V (2-10) 



2 

with the corresponding energies 

E sp in = ±^o-B? (2-H) 
to be added to the energy (|2.3j). Here, 



eh 

u = ^ 9.3 • 1(T 21 erg • Gauss" 1 = 5.8 ■ 10" 9 eV • Gauss" 1 (2.12) 

2mc 
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is Bohr magneton. Below, we do not consider the energy ( |2.11| ) related to the 
interaction of the spin of electron with the external magnetic field because 
it decouples from the orbital motion and can be accounted for separately. 
As to the numerical estimation, at B = B = 2.4 • 10 9 Gauss, the energy 
E spin ~ ±13.6 eV. 

For the ground level, i.e. at n = and s — 0, and zero momentum of the 
electron in the ^-direction, i.e. hk z = 0, we have from 



ehB 
2mc' 



(2.13) 



and due to Eq. (|2.6| ) the corresponding normalized ground state wave function 
is 

^ooo(r, <p, z) = ^ooo(r) = J^e^' 2 , (2.14) 



Io°Io n rdrdif |^ oo| 
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Figure 2: Landau ground state wave function of a single electron, ?/> oo (solid 
curve), Eq. ( 2.14 ), in strong external magnetic field B = Bq = 2.4 • 10 9 
Gauss, as a function of the distance r in cylindrical coordinates, and (for 
a comparison) the hydrogen ground state wave function (at zero external 
magnetic field), (1/ y/n)e~ r ^ a ° (dashed curve), as a function of the distance 
r in spherical coordinates. The associated probability densities are shown in 
Fig. |; 1 a.u. = a = 0.53 ■ 10" 8 cm. 
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Figure 3: The probability density for the case of the Landau ground state 
of a single electron, 27rr |-^ oo 1 2 (solid curve), Eq. ( |2.14j ), in the strong ex- 
ternal magnetic field B = B = 2.4 • 10 9 Gauss, as a function of the dis- 
tance r in cylindrical coordinates, and (for a comparison) probability den- 
sity of the hydrogen atom ground state (at zero external magnetic field), 
47rr 2 1(1/ \/7r)e~ r / ao | 2 (dashed curve), as a function of the distance r in spher- 
ical coordinates. The associated wave functions are shown in Fig. 1 a.u. 
= 0.53 • 10" 8 cm. 



The corresponding (smallest) Landau radius of the orbit of electron is 




(2.15) 



in terms of which ?/>ooo reads 

^»=S e ^- (2 - i6) 

Figure |] depicts ground state wave function of a single electron, ^ooo; 
in strong external magnetic field B = B = 2.4 • 10 9 Gauss (Rq = 1 a.u.), 
and (for a comparison) of the hydrogen ground state wave function, at zero 
external mag netic field, (1/ V?)e- r/a ° . Fi gures | and f| display the associated 
probability density of electron as a function of the distance r from the center 
of the orbit, the radius of which is about 1 a.u. 
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Figure 4: A contour plot of the (r, ip) probability density for the case of 
the Landau ground state of a single electron, 2nr\ip 000 | 2 , Eq. ( |2.14j ), in the 



strong external magnetic field B = B Q = 2.4 • 10 9 Gauss, as a function of the 
distance in a.u. (1 a.u. = 0.53 • 10~ 8 cm). Lighter area corresponds to higher 
probability to find electron. The set of maximal values of the probability 
density is referred to as an "orbit". 

The condition that the Landau radius Ro is smaller than the Bohr radius, 
i?o < clq, which is adopted here as the condition of a considerable "deforma- 
tion" of the electron orbit of the hydrogen atom, then implies 

2 '3 
Try CG 

B > B = — p— = 2.351 ■ 10 9 Gauss, (2.17) 

where m is mass of electron. Equivalently, this deformation condition corre- 
sponds to the case when the binding energy of the hydrogen atom, \EQ° hr \ = 
— me y2h 2 \ = 0.5 a.u. = 13.6 eV, is smaller than the ground Landau energy 

The above critical value of the magnetic field, Bq, is naturally taken as 
an atomic unit for the strength of the magnetic field, and corresponds to the 
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case when the pure Coulomb interaction energy of electron with nucleus is 
equal to the interaction energy of a single electron with the external magnetic 
field, \Eq oHt \ = Eq = 13.6 eV, or equivalently, when the Bohr radius is equal 
to the Landau radius, a = R = 0.53 • 10~ 8 cm. 

It should be stressed here that we take the characteristic parameters, 
Bohr energy \EQ° hr \ and Bohr radius do, of the hydrogen atom with a single 
purpose to establish criterium for the critical strength of the external mag- 
netic field, for the hydrogen atom under consideration. For other atoms the 
critical value of the magnetic field may be evidently different. 

3 Hydrogen atom in magnetic fields 

After we have considered quantum dynamics of a single electron in the ex- 
ternal magnetic field, we turn to consideration of the hydrogen atom in the 
external static uniform magnetic field. 

We choose cylindrical coordinate system (r,ip,z), in which the external 
magnetic field is B = (0,0,5), i.e., the magnetic field is directed along the 
z-axis. General Schrodinger equation for the electron moving around a fixed 
proton (Born-Oppenheimer approximation) in the presence of the external 
magnetic field is 

- 1. ^ + ;* + >i + * + ~ ^ + 2i ^) " = E ^ 

(3.1) 

where 7 = eH/2hc. Note the presence of the term 21^8^, which is linear in 
B, and of the term — 7 2 r 2 , which is quadratic in B. The other terms do not 
depend on the magnetic field B. 

The main problem in the nonrelativistic study of the hydrogen atom in 
the external magnetic field is to solve the above Schrodinger equation and 
find the energy spectrum. 

In this equation, we can not directly separate all the variables, r, if, and 
z, because of the presence of the Coulomb potential, e 2 / y/r 2 + z 2 , which does 
not allow us to make a direct separation in variables r and z. 

Also, for the general case of an arbitrary intensity of the magnetic field, 
we can not ignore the term 7 2 r 2 in Eq. (|3.1|) , which is quadratic in B, since 
this term is not small at high intensities of the magnetic field. 
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For example, at small intensities, B ~ 2.4 • 10 4 Gauss = 2.4 Tesla, the 
parameter 7 = eH/2%c ~ 1.7 • 10 11 cm~ 2 so that for (r) ~ 0.5 -10~ 8 cm we get 
the estimation 7 2 (r 2 ) ~ 10 6 cm" 2 and Landau energy hQ/2 is of the order 
of 10" 4 eV, while at high intensities, B ~ B = 2.4 • 10 9 Gauss = 2.4 • 10 5 
Tesla, the parameter 7 = 70 ~ 1.7 • 10 16 cm -2 so that the estimations are: 
7 2 (r 2 ) ~ 10 16 cm -2 (greater by ten orders) and the ground Landau energy is 
about 13.6 eV (greater by five orders). 

Below, we turn to approximate ground state solution of Eq. ( |3.1|) for the 
case of very strong magnetic field. 



3.1 Very strong magnetic field 

Let us consider the approximation of a very strong magnetic field, 

B > B = 2.4 • 10 9 Gauss. (3.2) 

Under the above condition, in the transverse plane the Coulomb interaction 
of the electron with the nucleus is not important in comparison with the 
interaction of the electron with the external magnetic field. So, in accord to 
the exact solution (|2.6Q for a single electron, one can seek for an approximate 



ground state solution of Eq. (|3.1| ) in the form (see, e.g. |J) of factorized 



transverse and longitudinal parts, 

if; = e-^/ 2 X (z), (3.3) 



where we have used Landau wave function (|2.14j) , s — 0, I — 0, and x( z ) 



denotes the longitudinal wave function to be found. This is so called adiabatic 
approximation. The charge distribution of the electron in the (r, ip) plane 
is thus characterized by the Landau wave function ~ e'^'\ i.e. by the 
azimuthal symmetry. 

In general, the adiabatic approximation corresponds to the case when the 
transverse motion of electron is totally determined by the intense magnetic 
field, which makes it "dance" at its cyclotron frequency. Specifically, the 
radius of the orbit is then much smaller than the Bohr radius, i?o ^ «0i or, 
equivalently, the Landau energy of electron is much bigger than the Bohr 
energy Eq ^> \EQ° hr \ = 13.6 eV. In other words, this approximation means 
that the interaction of electron with the nucleus in the transverse plane is 
ignorably small (it is estimated to make about 2% correction at B ~ 10 12 
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Gauss), and the energy spectrum in this plane is defined solely by the Lan- 
dau levels. The remaining problem is thus to find the longitudinal energy 
spectrum, in the z direction. 

Inserting wave function Q3.3|) into the Schrodinger equation (|3~T|), multi- 
plying it by if)*, and integrating over variables r and ip in cylindrical coordi- 
nate system, we get the following equation characterizing the z dependence 
of the wave function: 

(-I^ + $ +c m)* w ^ w ' (3 - 4) 

where 

/•oo p— P „ 

C(z) = -v^e 2 / dp = -e 2 ^^ erfc( v ^|z|), (3.5) 

JO y/p + 'JZ 2 

erfc(x) = 1— erf(x), and 




is error function. 

In Fig. [|, the exact potential C(z) is plotted, at B = 2.4 • 10 11 Gauss = 
100-Bo, for which case the lowest point of the potential is about C(0) — —337 
eV. In this case, in the (r, ip) plane, the hydrogen atom is characterized by 
the Landau energy Eq = 1360 eV ^> 13.6 eV, Landau radius R = 0.5 • 10~ 9 
cm <C 0.5 • 10 -8 cm, and the cyclotron frequency is Q = 4 • 10 18 sec -1 . For a 
comparison, the Coulomb potential —e 2 /\z\ is also shown. One can see that 
the Coulomb potential does not reproduce the effective potential C(z) to a 
good accuracy at small z. 

The arising effective potential C(z) is of a nontrivial form, which does 
not allow us to solve Eq. fl3.4p analytically. Below we approximate it by 
simple potentials, to make estimation on the ground state energy and wave 
function of the hydrogen atom. 

3.1.1 The Coulomb potential approximation 

At high intensity of the magnetic field, 7 ^> 1 so that under the condition 
7(z 2 ) > 1 we can ignore p in the square root in the integrand in Eq. ( |3.5| ). 
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Figure 5: The exact effective potential C(z) (top panel), at B — 2.4 • 10 11 
Gauss = 100-Bo, C(0) = —337 eV, and (for a comparison) Coulomb potential 
— e 2 /\z\ (bottom panel). 



Then, one can perform the simplified integral and obtain the result 



C(z) ~ V(z) = ~, at ^(z 2 ) > 1, (3.7) 
\z\ 

which appears to be a pure Coulomb interaction of electron with the nu- 
cleus, in the z direction. Therefore, Eq. ( |3.4| ) reduces to one- dimensional 
Schrodinger equation for Coulomb potential, 

;2 >2 „2 t2. 



h z d 2 e 2 ^ 7 



2m dz 2 



+ — + —L+E\ X (z)=0. (3.8) 



m 



13 



In atomic units (e = H = m = 1), using the notation 

E' = ^ + E, n 2 - 1 



m 



-2E n 



and introducing new variable 



x 



2z 



n 



we rewrite the above equation in the following form: 

d 



dx 2 







(4 









X(x) = 0, 



(3.9) 



(3.10) 



(3.11) 



where we assume, to simplify representation, x > 0. Introducing new func- 
tion v(x), 

x (x)=xe- x ' 2 v{x), (3.12) 
we get the final form of the equation, 



xv" + (2 — x)v ' — (1 — n)v = 0. 



(3.13) 



Note that this equation is identical to that obtained for the radial part of the 
three-dimensional Schrodinger equation for the hydrogen atom, at orbital 
quantum number zero. We note that it is a particular case of Cummer's 
equation, 

xv" + {b — x)v' — av = 0, 
the general solution of which is given by 



where 



and 



v(x) = C\ i-Fi (a, b, x) + C2U (a, b, x) 



MaAx) = v(h Vih \ v( . f\*H a -\l-t) 1 
T(b — a)T(a) Jo 



-a-l 



dt 



U(a,b,x) = — / e-*H a -\l + t) 1 
1 (a) Jo 



-a-l 



dt 



(3.14) 

(3.15) 
(3.16) 
(3.17) 



are the confluent hypergeometric functions, and are constants. In our 
case, a = 1 — n and b = 2. 
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For the wave function we then have 



X (x) = \x\e~^ /2 [Cf iFi(l - n, 2, \x\) + Of - n, 2, |x|)] , (3.18) 

where x given by Eq. ( |3.10| ) is now allowed to have any real value, and the " ±" 
sign in C^ 2 corresponds to the positive and negative values of x, respectively 
(the modulus sign is used for brevity). 
Now we turn to the normalization issue. 

The first hypergeometric function i-F\(l — n,2,x) is finite at x = for 
any n. At big x, it diverges exponentially, unless n is an integer number, 
n = 1,2,..., at which case it diverges polynomially. 

The second hypergeometric function U(l — n,2,x) behaves differently, 
somewhat as a mirror image of the first one. In the limit x — > 0, it is 
finite for integer n = 1,2,3, .. ., and diverges as 1/x for noninteger n > 1 
and for < n < 1. In the limit x — > oo, it diverges polynomially for 
integer n, tends to zero for noninteger n > 1 and for n = 0, and diverges for 
noninteger < n < 1. In Figs. [18] and [19] of Appendix, the hypergeometric 
functions ^(1 -n, 2, x) and £7(1 -n, 2, x) for n = 0,1/4,1/2,1,3/2,2,5/2,3 
are depicted (n may be taken, for example, n = 0.133 as well). 

In general, because of the prefactor xe~ x ^ 2 in the solution (|3.18|) which 



cancels some of the divergencies arising from the hypergeometric functions, 
we should take into account both of the two linearly independent solutions, 
to get the most general form of normalizable wave functions. 

As a consequence, the eigenvalues may differ from those corresponding 
to n = 1,2,... (which is an analogue of the principal quantum number in 
the ordinary hydrogen atom) so that n is allowed to take some non-integer 
values from to oo, provided that the wave function is normalizable. 

We are interested in the ground state solution, which is an even state. 
For even states, in accord to the symmetry of the wave function under the 
inversion z — > — z, we have 

Ct = C{, C+ = C^, x'(0) = 0. (3.19) 

Also, since n = 1 gives E' = — l/(2n 2 ) = —1/2 a.u., we should seek a 
normalizable wave function for n in the interval < n < 1, in order to 
achieve the energy value, which is lower than —1/2 a.u. If it is successful, 
the value n = 1 indeed does not characterize the ground state. Instead, it 
may correspond to some excited state. 
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One can see that the problem is in a remarkable difference from the ordi- 
nary three-dimensional problem of the hydrogen atom, for which the principal 
quantum number n must be integer to get normalizable wave functions, and 
the value n = 1 corresponds to the lowest energy. The distinction is due to 
the fact that the additional factor x in the solution ( |3.18| ) does not cancel in 
the total wave functionF]. 

We restrict consideration by the even state so we can focus on the x > 
domain, and drop the modulus sign, 

X (x) = xe~ x/2 [Cf - n, 2, x) + CfU{l - n, 2, x)] . (3.20) 

In summary, the problem is to identify a minimal value of n in the interval 
< n < 1, at which the wave function x( x ) represented by a combination of 
the two linearly independent solutions is normalizable and obeys the condi- 
tions ftnsp . 

The condition x'(0) = implies 

\e~ x / 2 [d ((2 - x^F^l - n, 2, x) + (1 - n)x 1 F 1 (2 -n,3,x)) + 

(3.21) 

+C 2 ((2 - x)U{l - n, 2, x) - 2(1 - n)xU{2 - n, 3, x))] lx=0 = 0. 

The function i-F\(l — n, 2, x) is well-defined at x = for any n. The function 
U(l — n,2,x) for small x behaves as a; _1 r _1 (l — n) so for x — it appears 
to be infinite, unless n = 1. The function U(2 — n, 3, x) for small x behaves 
as x" 2 r _1 (2 — n) so for x — it diverges (r(a) is Euler gamma function). 
Therefore, there is no way to satisfy the above condition unless we put C 2 = 
0, i.e. eliminate the second hypergeometric function, {7(1 — n,2,x), as an 
unphysical solution. As the result, only the first hypergeometric function 
determines the wave function, at which case it is well-defined only for integer 
n. 

For the first hypergeometric function as the remaining solution, the situ- 
ation is well known, 

iFi(a, 6, 1x1)^=0 = 1, ^ 1 F x {a,b,x)\ x=Q = ^. (3.22) 



the three-dimensional Coulomb problem it does cancel in the total wave function. 
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We take n = 1 as this value implies the lowest energy E' in this case. The 
energy for n = 1 is E' = — 1/2 a.u., i.e. 

2ft 2 2 v ; 



where we have used h j/m = due to Eqs. Q2.2p and 

Thus, the particular state n = 1 of the hydrogen atom is characterized by 
usual Bohr distance ao in the z direction, with the longitudinal wave function 

X(z) ~ \z\e- lzl/ao . (3.24) 

In addition, it is easy to check directly, without using the general hypergeo- 
metric function technique, that x( z ) — ze ~ z satisfies the equation (^d 2 /dz 2 + 
l/z + e)x{z) = 0, for e = -1/2. 

In the (r, ip) plane (n = 0, s = 0) the orbit is characterized by the Landau 
radius Ro -C ao, with the Landau wave function ~ e~ r l 4R a, and the Landau 
energy Eq given by the second term in Eq. ( |3.23| ). Thus, the n = 1 wave 
function in the Coulomb approximation is written as 



2 ^-7r 2 /2-kl/ao 



il)(r,<p,z) ~ y^e" 7r ^-l-IM) = 4ii o a o . (3.25) 

Also, spin of the electron for the ground state is aligned antiparallel to the 
magnetic field. 

The associated probability density is evidently of a cylindrical (axial) 
symmetry and can be described as two Landau orbits of radius Rq in different 
(r, ip) planes, one at the level z = —a Q , and the other at the level z = +a Q , 
with the nucleus at z — 0, as schematically depicted in Fig. EL Presence of 



two Landau orbits occurs in accord to the wave function ( |3.25| ) , which equals 
zero at z = and is symmetrical with respect to the inversion, z — > —z. The 
electron moves simultaneously on these two orbits. 

The above electron charge distribution supports the study made by San- 
tilli who proposed the polarized toroidal electron orbit in the hydrogen 
atom under the action of strong magnetic field, but the Coulomb approxi- 
mation suggests the electron charge distribution in the form of two identical 
coaxial round orbits separated by relatively big vertical distance; see Fig. |. 

Note that the size of the hydrogen atom in the z direction is predicted 
to be about 10 times bigger than that in the transverse plane, i.e. Rq = 
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Figure 6: Schematic view on the hydrogen atom in the ground state, at very 
strong external magnetic field B = (0, 0, B), B S> B = 2.4 • 10 9 Gauss, due 
to the Coulomb approximation approach, Eq. ( |3.25| ). One electron moves 
simultaneously on two Landau orbits of radius Ro which are shown schemat- 
ically as torii in the different (x, y) planes, one torus at the level z = — 
and the other at the level z = +ao, with the nucleus shown in the center at 
z — 0. Each torus represents the (x, y) probability distribution as shown in 
Fig. H but with small Landau radius, Ro <C a . The spin of the electron is 
aligned antiparallel to the magnetic field. 

0.53 ■ 10~ 9 cm ~ O.lao; the hydrogen atom is thus highly elongated in the z 
direction (this ratio is not kept in Fig. The cyclotron frequency ( |2.4j ) is 
of about 4 • 10 18 sec" 1 . Classically, the linear velocity of the electron on such 
an orbit is v = 2ttRqQ and is estimated to be v = 1.3 • 10 10 cm/sec ~ 0.4c. 

We should to note here that Landau radius R , at B = 2.4 • 10 11 Gauss 
= 100i?o, is about ten times bigger than Compton wavelength of electron, 
A e /27r ~ 0.4- 10~ 10 cm, hence we can ignore quantum electrodynamics effects. 

Finally, it should be noted that the above used Coulomb approximation 
is rather a crude one (see Fig. [5]). Indeed, the energy in the ^-direction does 
not depend on the magnetic field intensity while it is obvious that it should 
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depend on it. In the next Section, we approximate the effective potential 
C(z) with a better accuracy. 



3.1.2 The modified Coulomb potential approximation 

We note that due to the exact result Q3.5|), the effective potential C(z) tends 
to zero as z — > oo (not a surprise). However, a remarkable implication of the 
exact result is that C(z) is finite at z — 0, namely, 

C(0) = -v^7e 2 , (3.26) 



so that the effective potential C(z) indeed can not be well approximated 
by the Coulomb potential, —e 2 /\z\, at small z (\z\ < 1CT 8 cm) despite the 
fact that at big values of 7 (7 ^> 70 — 1.7 ■ 10 16 cm -2 ) the Coulomb po- 
tential reproduces C(z) to a good accuracy at long and medium distances. 



Thus, the approximation considered in Sect. |3.1.1| appears to be not valid at 
(important) short distances. 

Consequently, Eq. ( |3.4| ) with the exactly calculated C(z) will yield the 
ground state of electron in the z direction which is drastically different from 
that with the Coulomb potential considered in Sect. |3.1.1| . Thus, very intense 
magnetic field essentially affects dynamics of electron not only in the (r, tp) 
plane but also in the z direction. We expect that for the case of exact 
potential (|3.5|) the ground state is characterized by much lower energy in the 
z direction as compared to that of the one-dimensional Bohr state ( |3.23| ). 

The exact potential C(z) can be well approximated by the modihed 
Coulomb potential, 

C(z) ~ V{z) = -t-£— , (3.27) 

\Z\ + Zo 

where z Q is a parameter, z 7^ 0, which depends on the field intensity B due 
to 

e 2 1 / 2hc 



The analytic advantage of this approximation is that V(z) is finite at z = 0, 
being of Coulomb-type form. 

Calculations for this approximate potential can be made in an essentially 
the same way as in the preceding Section, with \x\ being replaced by \x\ +Xq, 



19 



where xq > (we remind that x = 2z/n), 

X (x) = (|x|+x )e- (|:c|+:co)/2 [CtxFi(l-n,2, \x\+x ) +C£U(l-n,2, \x\+x ) 

(3.29) 

An essential difference from the pure Coulomb potential case is that the 
potential ( p.27| ) has no singularities. The x — > asymptotic of the associated 
second hypergeometric function, U(l — n, 2, \x\ + Xo) becomes well-defined for 
noninteger n. Thus, the condition x'( x )\x=o — does not lead to divergencies 
in the second solution for noninteger n < 1. It remains only to provide well- 
defined behavior of the wave function at x — > oo. 

Analysis shows that normalizable wave functions, as a combination of two 
linearly independent solutions, for the modified Coulomb potential does exist 
for various non-integeif\ n. We are interested in the ground state solution, 
so we consider the values of n from to 1. Remind that E' = —l/(2n) 2 so 
that for n < 1 we obtain the energy lower than —0.5 a.u. 

For n < 1, the first hypergeometric function is not suppressed by the 
prefactor xe~ x l 2 in the solution (|3.29 ) at large x so we are led to discard it as 



an unphysical solution by putting C\ — 0. A normalizable ground state wave 
function for n < 1 is thus may be given by the second term in the solution 
( |3.29|) . The condition x'(x)\x=o = implies 



l e -( x + x °)/ 2 C 2 [(2 -x- x )U(l -n,2,x + x )- 

(3.30) 

-2(1 - n)(x + x )U(2 — n,3,x + x ))]\ x=0 = 0. 

The l.h.s of this equation depends on n and x , so we can select some field 
intensity B, calculate associated value of the parameter xq = Xq(B) and find 
the value of n, from which we obtain the ground state energy E' . 

On the other hand, for the ground state this condition can be viewed, 
vice versa, as an equation to find x for some given value of n |J. Hereby we 
reverse the order of the derivation, to simplify numerical calculations. 

For example, taking the noninteger value 



n = 1/V15.58 ~ 0.253 < 1 (3.31) 
we find from Eq. ( p.30|) numerically 

x = 0.140841. (3.32) 

2 This is the case when one obtains a kind of "fractional" quantum number n which 
does make sense. 



20 



This value is in confirmation with the result xq = 0.141 obtained by Heyl and 
Hernquist (see Discussion). On the other hand, xq is related in accord to 
Eq. ( |3.28| ) to the intensity of the magnetic field, 

2z I 8hc 
%o = — = \ — 5-^, 3.33 

from which we obtain B ~ 4.7 • 10 12 Gauss. Hence, at this field intensity 
the ground state energy of the hydrogen atom is given by E' = — l/(2n 2 ) = 
— 15.58 Rydberg. 

n=0.253, x0=0.141 

0.8 
0.6 
o 0.4 
0.2 


-10 -5 5 10 

x 




Figure 7: The longitudinal ground state wave function Eq. (|3.34|) , of 

the hydrogen atom in the mag netic field B = 4.7 • 10 12 Gauss; n = l/y/^2E, 
x = 2z/n. 

The longitudinal ground state wave function is given by 

X(x) ~ (|x| + x )e (|:E|+a;o)/n £/(l - n, 2, |x| + x ), (3.34) 
and is plotted in Fig. 0). The total wave function is 



x {x) ~ J— m e "%{\x\ + x )e^ +xo)/n U(l - n, 2, |x| + x ), (3.35) 
y ItiKq 

and the associated three-dimensional probability density is schematically de- 
picted in Fig. ||. 

In contrast to the double Landau-type orbit implied by the Coulomb 
potential approximation, the modified Coulomb potential approach provides 
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Figure 8: Schematic view of the hydrogen atom in the ground state, at very 
strong external magnetic field B = (0,0, B), B ^> B = 2.4 ■ 10 9 Gauss, 
due to the modified Coulomb approximation approach. The electron moves 
on the Landau orbit of small radius i?o ^ 0.53 ■ 10 -8 cm, which is shown 
schematically as a torus. The vertical size of the atom is comparable to i?o- 
Spin of the electron is aligned antiparallel to the magnetic field. 

qualitatively correct behavior and much better accuracy, and suggests a single 
Landau-type orbit shown in Fig. |8] for the ground state charge distribution 
of the hydrogen atom. This is in full agreement with Santilli's study || of 
the hydrogen atom in strong magnetic field. 

Only excited states are characterized by the double Landau-type orbits 
depicted in Fig. [| We shall not present calculations for the excited states in 
the present paper, and refer the reader to Discussion for some details on the 
results of the study made by Heyl and Hernquist ||. 

3.1.3 Variational and numerical solutions 

Review of approximate, variational, and numerical solutions can be found 
in the paper by Lai [|J. Accuracy of numerical solutions is about 3%, for 
the external magnetic field in the range from 10 11 to 10 15 Gauss. Figure ||] 
reproduces Fig. 1 of Ref. M. It represents the computed ionization energy 
Qi of the hydrogen atom and dissociation energy of H 2 molecule as functions 
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Figure 9: Ionization energies of the hydrogen atom and H 2 molecule exposed 
to very strong magnetic field; B u = B ■ 10~ 12 Gauss (reproduction of Fig. 1 
by Lai §). 



of the intensity of the magnetic field; B u = B ■ 10 12 Gauss, i.e. Bi 2 = 1 
means B = 10 12 Gauss. 




Figure 10: Longitudinal z-size of the hydrogen atom in the magnetic field, Eq. 
( |3.36| ), as a function of the magnetic field intensity; b = B / B , B = 2.4 • 10 9 
Gauss, 1 a.u. = 0.53 ■ 10~ 8 cm. 

Due to the variational results || , the z-size of the hydrogen atom in the 
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Figure 11: Transverse size of the hydrogen atom in the magnetic field, Eq. 
( |3.37| ), as a function of the magnetic field intensity; b = B / B , B = 2.4 • 10 9 
Gauss, 1 a.u. = 0.53 ■ 10~ 8 cm. 



ground state is well approximated by the formula (Fig. pTDf ) 

1 



MB /B ) 

the transverse (Landau) size is (Fig. [11]) 

1 



-a 



a.u.; (3.36) 



a.u.; (3.37) 



'B/B 

and the ground state energy is (Fig. [12]) 

E ~ -0.16[MB/B )] 2 a.u., (3.38) 

with the accuracy of few percents, for b = B / B in the range from 10 2 to 
10 6 . 

One can see for B = 100B , that the variational study predicts the ground 
state energy E = —3.4 a.u. = —92.5 eV (much lower than the value, —0.5 
a.u., predicted by the n = 1 Coulomb approximation in Sect. |3.1.1| ), the 
transverse size L± of about 0.1 a.u. (the same as in the Coulomb approxi- 
mation) and the z-size L z of about 0.22 a.u. (much smaller than the value, 
1 a.u., predicted by the Coulomb approximation). This confirms the result 
of the approximate analytic approach, in which n < 1 implies ground state 
energy much lower than —0.5 a.u. 
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Figure 12: Variational ground state energy of the hydrogen atom in the 
magnetic field, Eq. ( |3.38| ), as a function of the magnetic field intensity; 
b = B/B , B = 2.4 ■ 10 9 Gauss, 1 a.u. = 27.2 eV. 

4 Discussion 

In a physical context, the used adiabatic approximation implies that the 
position of the nucleus is not "fixed" in the (r,ip) plane. Indeed, there is no 
reason for the nucleus to stay "exactly" at the center of the orbit when the 
Coulomb force in the transverse plane is totally ignored. 

In contrast, in the z direction, we have more complicated situation, 
namely, the magnetic field leads to the effective interaction potential C(z), 
which is finite at z = and is Coulomb-like in the long-range asymptotic. 
Hence, the position of the nucleus remains to be fixed in the z direction 
relative to the orbit of the electron (bound state). 

Due to the (crude) Coulomb approximation of Sect. |3.1.1| , the electron 
probability density of the hydrogen atom under very intense magnetic field, 
at n — 1, appears to be of rather nontrivial form schematically depicted in 
Fig.|. 

More accurate approach may yield some different electron charge distri- 
butions but we expect that the symmetry requirements (namely, the axial 
symmetry and the (anti) symmetry under inversion z — » —z) would lead again 
to a toroidal orbit, which is splitted into two identical coaxial orbits (the torii 
are separated) similar to that shown in Fig. ^|, or to an unsplitted one (the 
two torii are superimposed) as shown in Fig. || 

Accurate analytic calculation of the ground and excited hydrogen wave 
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functions made by Heyl and Hernquist || in the adiabatic approximation 
leads to the longitudinal parts of the wave functions (the vertical z depen- 



dence) shown in Fig. O, which reproduces the original Fig. 3 of their work; 
( = 27iaz/X e ; B = 4.7- 10 12 Gauss. They used the modified Coulomb poten- 
tial of the type (|3.27|) , and the additional set of linearly independent solutions 



(|x| + x m )e-^l +a; -)/ 2 1 F 1 (l - n, 2, |x| + x m ) f 



■ dt. 



of the one-dimensional modified Coulomb problem in the form 

(t 1 F 1 (l-n,2,t)y 

(4.1) 

where m = corresponds to the ground state. For the ground state with n = 
1/V15.58, i.e. the binding energy is 15.58 Rydberg, they found x = 0.141, 
which corresponds to B = 4.7 • 10 12 Gauss. This result is in agreement with 
the study made in Sect. |3.1.2j . 




Figure 13: The axial wavefunctions of hydrogen in an intense magnetic field 
(analytic calculation) for B = 4.7- 10 12 Gauss. The first four even states with 
axial excitations, 1 000) (ground state), 1 002) , 1 004) , and 1 006) (left panel), and 
odd states 1 001) and 1 003) (right panel) are depicted. Here, B = 4.7 • 10 12 
Gauss, ground state energy is E = -15.58 Rydberg, n = l/y/l5l8, ( = 
2z/n corresponds to x in our notation; z in a.u., 1 a.u. = 0.53 • 10~ 8 cm 
(reproduction of Figure 3 by Heyl and Hernquist ||). 
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Figure 14: The axial ground state 1 000) wavefunction of hydrogen in an 
intense magnetic field, analytic calculation (solid curve) and numerical result 
(dashed line); B = 4.7 • 10 12 Gauss, E = -15.58 Rydberg, n = 1/Vl^58, 
( = 2z/n corresponds to x in our notation; z in a.u., 1 a.u. = 0.53 • 10~ 8 cm 
(reproduction of Figure 4 by Heyl and Hernquist ||). 

One can see from Fig. [13] that the peak of the ground state wave function 
1 000) is at the point z = (see also Fig. [L4| for a close view) while the largest 
peaks of the excited wave functions are away from the point z = (as it 
was expected to be). Consequently, the associated longitudinal probability 
distributions (square modules of the wave functions multiplied by the volume 
factor of the chosen coordinate system) are symmetric with respect to z — > 
—z, and their maxima are placed in the center z = for the ground state, 
and away from the center for the excited states. 

In contrast to the double toroidal orbit (Fig. ^[) implied by the Coulomb 
approximation of Sect. [3.1.1| , the modified Coulomb potential approach sug- 
gests a single toroidal orbit (Fig. H) for the ground state charge distribution 
of the hydrogen atom. This confirms Santilli's study @ of the hydrogen 
atom in strong magnetic field. The excited states \00u) are characterized by 
a double toroidal orbit schematically depicted in Fig. 

The computed ground state 1 000) binding energy of the hydrogen atom 
for different field intensities are ||: 
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Magnetic field B 


Binding energy, OOO) state 


(Gauss) 


(Rydberg) 


4.7 x 10 12 


15.58 


9.4 x 10 12 


18.80 


23.5 x 10 12 


23.81 


4.7 x 10 13 


28.22 


9.4 x 10 13 


33.21 


23.5 x 10 13 


40.75 


4.7 x 10 14 


47.20 



They calculated first-order perturbative corrections to the above energies 
and obtained the values, which are in a good agreement with the results by 
Lai ]§ (see Fig. |) and Ruder et al. 0. 

It is remarkable to note that in very strong magnetic field, the energy 
differences between the excited states and the ground state of the hydrogen 
atom are small in comparison to the absolute value of the ground state energy. 
The excited states are represented by the double torii configurations since 
the peaks of the associated wave functions spread far away from the center 
z = in both +z and — z directions. 

The above ground state electron charge distributions strongly support the 
study made by Santilli who proposed the polarized toroidal electron orbit 
in the hydrogen atom under the action of strong magnetic field as a physical 
picture at the foundation of the new chemical species of magnecules. The 
only addition to this picture is that the the excited states are characterized 
by the double torii configuration. 

Since zero field ground state case is characterized by perfect spherically 
symmetric electron charge distribution of the hydrogen atom, intermediate 
intensities of the magnetic field are naturally expected to imply a distorted 
spherical distribution. However, deeper analysis is required for the interme- 
diate magnetic field intensities because the adiabatic approximation is not 
longer valid in this case. Namely, an interesting problem is to solve the 
Schrodinger equation (|3.1|) for the case when intensity of the magnetic field 
is not very strong, B < Bq, with the corresponding characteristic Landau 
ground state energy of about the Bohr energy. In this case, both the Coulomb 
and magnetic interactions should be taken at the same footing that leads to 
complications in its analytic study. 
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Figure 15: Ionization energies of He atom exposed to very strong magnetic 
field (numerical solutions); B 12 = B ■ 10~ 12 Gauss (reproduction of Fig. 5 by 
Lai i). 

As to the multi-electron atoms, an interesting problem is to study action 
of a very strong external magnetic field on He atom (see. e.g., Heyl and Hern- 
quist || and Lai ||; Fig. [L5D and on the multi-electron heavy atoms, with 
outer electrons characterized by a nonspherical charge distribution, such as 
the p-electrons in Carbon atom, orbitals of which penetrate the orbitals of 
inner electrons. Very intense magnetic field would force such outer electrons 
to follow small round Landau orbits. In addition to the effect of a direct 
action of the magnetic field on the inner electrons, a series of essential re- 
arrangements of the whole electron structure of the atom seems to be occur 
with a variation of the external magnetic field strength. Indeed, the mag- 
netic field competes with both the Coulomb energy, which is different for 
different states of electrons, and the electron-electron interactions, including 
spin pairings. However, it is evident that at sufficiently strong fields, all the 
electron spins are aligned antiparallel to the magnetic field — fully spin po- 
larized configuration — while at lower field intensities various partial spin 
polarized configurations are possible. 
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Figure 16: Contour plots of the (r, z) plane electron density of the iron atom 
according to density matrix theory at two different magnetic field strengths, 
10 11 Gauss (left) and 10 12 Gauss (right). The outermost contour encloses 99% 
of the negative charge, the next 90%, then 80% etc., and the two innermost 
5% and 1% respectively (reproduction of Fig. 5 by Johnsen and Yngvason 

§)■ 

In accord to the numerical calculations based on the density matrix the- 
ory approach by Johnsen and Yngvason |J, which is in quite good agreement 
with the result of Hartree-Fock approach at very strong magnetic field inten- 
sities, the inner domain in the iron atom (26 electrons) is characterized by 
slightly distorted spherically symmetric distribution even at the intensities 
as high as B = 100£>o • • • 1000i?o- The outer domain appears to be of spe- 
cific, highly elongated distribution along the direction of the magnetic field 
as shown in Fig. [16]. The possible interpretation that the inner electrons re- 
main to have spherical distribution while outer electrons undergo the squeeze 
seems to be not correct unless the spin state of the iron atom is verified to 
be partially polarized. So, we can conclude that all the electrons are in 
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highly magnetically polarized state (Landau state mixed a little by Coulomb 
interaction), and the electron structure is a kind of Landau multi-electron 
cylindrical shell, with spins of all the electrons being aligned antiparallel to 
the magnetic field (fully spin polarized configuration). 

Another remark regarding Fig. [16] is that the contours indicating nearly 
spherical distribution will always appear since the Coulomb center (nucleus) 
is not totally eliminated from the consideration (non-adiabatic approxima- 
tion), and it forces a spherical distribution to some degree, which evidently 
depends on the distance from the center (closer to the center, more spheric- 
ity). We note that Fig. |1^ is in qualitative agreement with Fig. |6] in the sense 
that the predicted charge distribution reveals symmetry under the inversion 
z — > —z, with the characteristic ^-elongated toroidal orbits. 

An exciting problem is to study H 2 molecule under the action of strong 
external static uniform magnetic field using Schrodinger equation. Since the 
binding energy of H 2 molecule is of about 4.75 eV (much less than 13.6 eV of 
individual hydrogen atoms) and two interacting electrons of spin 1/2 partic- 
ipate dynamics we expect interesting physical predictions. For example, it 
is interesting to analyze the problem of bonding, in the presence of a strong 
magnetic field. 

However, before studying H 2 molecule it would be very useful to investi- 
gate much simpler two-center system, ff 2 ion, as it can give valuable infor- 
mation on the features of a two-center system under the action of a strong 
magnetic field. We refer the interested reader to Refs. @, ^] for details of 
studies on ion and H 2 molecule in strong magnetic field. Figure ^ dis- 
plays ground and first excited state wave functions of h 7 ^ i° n calculated by 
Heyl and Hernquist 0. 

Finally, it should be noted that studies on the atomic and molecular 
systems in strong magnetic fields and related issues are becoming extensive 
last years. We refer the reader to a recent review by Lai 0] for an extensive 
list of references, and below present summary of some selected recent papers. 

The work by Kravchenko et al. |J uses an orthogonal approach and 
achieves very high precision by solving the general problem of the hydrogen 
atom in an arbitrarily strong magnetic field. This analytical work is of much 
importance as it provides solid ground to analyze few-electron and multi- 
electron atoms, and simple diatomic molecules. 

A series of papers by Schmelcher et al. @ is devoted to a systematic study 
of various simple atomic, ionic and molecular systems in strong magnetic 
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Figure 17: The ground and first-excited state of H2" ion. The solid line traces 
1 000) , and the dashed line follows |0-10). The triangles give the positions 
of the protons for the ground state and the squares for the excited state. 
Magnetic field B = 4.7 • 10 12 Gauss is pointed along the internuclear axis; 
( = 2naz/X e denotes z in a.u.; 1 a.u. = 0.53 • 10~ 8 cm (reproduction of 
Figure 5 by Heyl and Hernquist ||). 



fields via Hartree-Fock method. One of the papers by Schmelcher et al. 
considers the carbon atom. 

Analytical approximations have been constructed by Potekhin for 



binding energies, quantum mechanical sizes and oscillator strengths of main 
radiative transitions of hydrogen atoms arbitrarily moving in magnetic fields 
~ 10 12 . . . 10 13 Gauss. This approach is of much importance in investigating 
an interaction of highly magnetized hydrogen atoms in a gaseous phase. 

Computational finite element method for time-dependent Schrodinger 
equation has been developed by Watanabe and Tsukada []TTJ to study dynam- 
ics of electrons in a magnetic field that can be used effectively for mesoscopic 
systems, such as small clusters. 
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Appendix 




Figure 18: The confluent hypergeometric functions iFi(l — n, 2, x) and U(l — 
n, 2, x) for various n. 



34 




Figure 19: The confluent hypergeometric functions iFi(l — n, 2, x) and U(l — 
n, 2, x) for various n. 
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